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ABSTRACT 

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts 
(dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this 
measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in 
very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower 
at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys 
over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear 
evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination 
of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated. 

Key words. Surveys - - Radio continuum: general - - Galaxies: active 



1. Introduction 

Planck (jTauber et al.ll2010t iPlanck Collaboration! [20TTah 
is the third-generation space mission to measure the 

1 Planck ( |http: / /www .esa.int /Planck ) is a project of the 
European Space Agency (ESA) with instruments provided by 
two scientific consortia funded by ESA member states (in par- 
ticular the lead countries: France and Italy) with contributions 
from NASA (USA), and telescope reflectors provided in a collab- 
oration between ESA and a scientific consortium led and funded 
by Denmark. 
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anisotropy of the cosmic microwave background (CMB). 
It observes the sky in nine frequency bands cov- 
ering 30-857 GHz with high sensitivity and angu- 
lar resolution fr om 31' to 5'. The Low Frequency 
Instr u ment (LFI; Mandolesi et all 120101 : iBersanelli et al.l 
I2010t iMennella et al.l 1201 lft covers the 30, 44, and 
70 GHz bands with amplifiers c ooled to 20 K . The 
High Frequency Instrument (HFI; iLamarre" etaLl [20101 
IPlanck HFI Core Team! |2011al) covers the 100, 143, 217, 
353, 545, and 857 GHz bands with bolometers cooled to 
0.1K. Polarization is measured in all but th e highest two 
bands (|Leahv et al.ll2010t iRosset et alj |2010). A combina- 
tion of radiative cooling and three mechanical coolers pro- 
duces the temperatures neede d for the detectors and optics 
(|Planck Collaborationll20lfbl ) . Two data processing centres 
(DP Cs) check and calibrate the data and make maps of th e 
sky (IPlanck HFI Core Team! 120111* IZacchei et &U [20ll . 
Planck's sensitivity, angular resolution, and frequency cov- 
erage make it a powerful instrument for galactic and ex- 
tragalactic astrophysics as well as cosmology. Early astro- 
physics results are given in Planck Collaboration, 2011h-z. 

The P lanck Early Release Compact Source Catalogue 
(ERCSC, IPlanck Collaboration! (|2011dl ) reports data on 
sources detected during the first 1.6 full-sky surveys, and 
thus offers, among other things, the opportunity of study- 
ing the statistical properties of extragalactic sources over 
a broad frequency range never fully explored by blind sur- 
veys. We will focus here on counts of extragalactic radio 
sources and on their spectral properties in the 30-217 GHz 
range Q 

Although knowledge of the statistical properties at 
high radio frequency for this population of extragalactic 
sources has greatly improved in the recent past - thanks 
to many ground-based observational campaigns and to the 
Wilkinson Microwave Anisotropy Probe ( WMAP) surveys 
from space - above about 70 GHz these properties are still 
largely unknown or very uncertain. This is essentially due 
to the fact that very large area surveys at mm wavelengths 
are made difficult by the small fields of view of ground- 
based radio telescopes and by the long integration times 
required. 

The most recent estimates on source number counts up 
to ~ 50 — 70 GHz, and the optical identifica tions of the cor- 
responding b right point sources (see, e.g., iMassardi et al.l 
(2008, 2010)), show that these counts are dominated by 
radio sources whose average spectral index is "flat", i.e., 
a ~ 0.0 (with the usual convention S v oc v a ). This re- 
sult confirms that the underlying source population is es- 
sentially made of Flat Spectrum Radio Quasars (FSRQ) 
and BL Lac objects, collectively called blazarsQ with mi- 
nor contributions coming from other sou rce populations 
(|Toffolatti et alJll998t Ide Zotti et al.ll2005r i. At frequencies 



2 In all our calculations we have us ed the effective cen- 
tral frequencies for the Planc k channels (|Mennella et al.ll2011l : 
IPlanck HFI Core Te am 2011a), although we indicate their nom- 
inal values. The most relevant difference is at 30 GHz, with a 
central frequency of 28.5. 

3 Blazars are jet-dominated extragalactic objects character- 
ized by a strongly variable and polarized emission of the non- 
thermal radiation, from low radio energi es up to the high en- 
ergy gamma rays (|Urrv fc Padovani 1995|). Detailed analyses of 
Spectral Energy Distributions (SEDs) of complete blazar sam- 
ples built by using simultaneous Planck , Swift and Fermi data 
are given in l|Planck C olfabora tionll2011lJ) . 



> 100 GHz, however, there is now new information for 
sources with flux densities below about 1 Jy coming from 
the S outh Pole Telescope (SPT) collaboration ()Vieira et al.l 
l2010h. with surveys over 87 deg 2 at 150 and 220 GHz, and 
from the Atacama Co smology Telescope (A CT) survey over 
455 deg 2 at 148 GHz (iMarriaee et al.ll2010h . 

The "flat" spectra of blazars are generally be- 
lieved to result from the superposition of different com- 
ponents in the inner part of AGN relativistic jets, 
each w ith a different synchrotron self -absorption fre- 
quency (|Kellermann fc Paulinv-Tothlll969D . At a given fre- 
quency, the observed flux density is thus dominated by 
the synchrotron-emitting component which becomes self- 
absorbed and, in the equipartition regime, the resulting 
spectrum is approximately flat. However, this "flat" spec- 
trum cannot be maintained up to very high frequencies, be- 
cause of electron energy losses in the dominant jet-emission 
component (i.e., electron ageing), or the transition to the 
optically-thin regime, with the onset of a "steep" spec- 
trum with a standard spectral index a = —0.7 to —0.8. 
A slightly steepened spectrum may also be caused by the 
superposition of many jet components. The redshift moves 
the observed steepening to lower frequencies and, thus, a 
greater fraction of blazar sources are observed with a steep 
spectrum at sub-mm wavelengths. With current data it 
is not yet possible to decide among the different scenar- 
ios. However, given their sensitivity and full sky coverage, 
Planck surveys are uniquely able to shed light on this tran- 
sition from an almost "flat" to a "steep" regime in the spec- 
tra of blazar sources. 

The outline of this paper is as follows. In i) l2.1l we briefly 
sketch the main properties of the ERCSC. In ^ 12.21 we sum- 
marize the source validation. In § 12.31 we describe the com- 
plete sample selected at 30 GHz, used for the analysis of 
spectral properties. In §[3]we present the source counts over 
the frequency range 30-217 GHz. In §U]we investigate the 
spectral index distributions in different frequency intervals. 
Finally, in §[5] we summarize our main conclusions. 



2. The Planck ERCSC 
2.1. Overview 



The Planck ERCSC (jPlanck Collaborationll2011cD lists po- 
sitions and flux densities for the compact sources recovered 
from the Planck first 1.6 full sky survey maps in nine fre- 
quency channels between 30 and 857 GHz. Thus about 60% 
of the sources have been covered twice, with a time sepa- 
ration of about 6 months. Sources near the ecliptic poles, 
where the scan circles intersect, are often covered multiple 
times. ERCSC flux densities are therefore averages over dif- 
ferent observing time periods. They have been calculated 
by aperture photometry, usi ng the most recent definition of 
the beam shapes and sizes (Pla nck HFI Core Teanj|2011bt 
IZacchei et al.ll201~il) . In the frequency range considered in 
this paper, the Planck photometric calibration is based on 
the CMB dipole and on the modulation ind uced on it by 
the sp acecraft orbital motion. According to IZacchei et al.l 
(|201lf ) the absolute photometric calibration at LFI fre- 
quencies (30, 44 and 70 GHz) is at the ~ 1% level, while 
IPlanck HFI Core Team! (|2011bl ) reports on a relative photo- 
metric accuracy, between the frequency channels from 100 
to 353 GHz, better than 2% and, more likely, at the ~ 1% 
level. 
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Fig. 1. Comparison between the ERCSC flux densities at 30 GHz (left panel) and at 44 GHz (right panel) with the 
almost simultaneous ATCA measurements (PACO project) at 32.2 and 39.7 GHz, respectively. No correction for the 
slightly different frequencies has been applied. 



The final versi o n of the ERCSC (see also 
iPlanck Collaboration! l2011d for more details) has 
been created by specifically applying the so-ca lled 
"PowellSnakes" detection method (|Carvalho et al.l[2009T ) to 
the Planck full sky (dipole subtracted) anisotropy maps 
in the frequency ch annels from 30 to 143 G Hz and the 
SExtractor package (jBertin fc Arnouti Il996f) in the four 
channels at higher frequencies (from 217 to 857 GHz). A 
"compact" source is accepted in the ERCSC if it survived 
a set of primary and secondary selection criteria. The 
primary criterion utilized the feedback from the Monte 
Carlo Quality Assessment system and introduced a signal- 
to-noise ratio cut to ensure that > 90% of the sources in 
the catalogue are reliable and have a flux density accuracy 
better than 30%. The secondary criterion comprised a 
set of cuts that removed the extended sources with an 
elongation > 3.0 (the ratio between the major and minor 
axis, in pixels, of the detected source) and sources that 
could be potentially spurious. For instance a source is 
dropped whenever > 5% of the pixels within 2xFWHM of 
its position had invalid values (see the ERCSC Explanatory 
Supplement for more details). The overriding requirement 
in constructing the ERCSC is source reliability and not 
completeness. 



2.2. Validation and photometry check of ERCSC sources 

The validation process for the compact sources included 
in the Planck ERCSC was performed by two different 
teams, selected among members of the Planck Consortia: a 
Validation Team on radio sources (VTRS) and a Validation 
Team on far-IR sources. The two teams worked separately 
at first, but cross-checked their results in the second phase 
of the process. The processing steps and main outcomes 
are summarized in the Explanatory Supplement released 
with the ERCSC. For compact ra dio (i.e., synchrotron 
domin ated) sources, the VTRS (see IPlanck Collaboration! 
(I2011ih for a more detailed discussion) has found that 
> 97% of the ERCSC sources at 30 GHz have reli- 
able counterparts in published catalogues at GHz fre- 



quencies CPM N: IWright et all 119961: GB6: iGreeorv et a! 
[19961: NVSS: ICondon et al.l 119981: SUMSS: iMauch et al. 
booa AT20G: iMassardi et all 120081: iMurphv et all l2010h . 
Similar (although slightly lower) percentages were found 
for ERCSC sources detected at 44 and 70 GHz. At higher 
frequencies (> 100 GHz) Planck detects an increasing frac- 
tion of dusty galaxies, undetected by low-frequency sur- 
veys. Therefore, the source reliability was confirmed by 
internal matches of sources detected in two neighboring 
Planck frequency channels: i.e., 143 and 217 GHz, or 217 
and 353 GHz, etc.. However, the validation of synchrotron- 
dominated sources is still relatively easy to perform, since 
all of them must be present in low- frequency catalog ues. 

The WMAP 7-year catalogue (jGold et all l201oj) con- 
tains a total of 471 sources detected in at least one fre- 
quency channel. Of these, 289, 281, 166 and 59 sources 
are detected as > 5a peaks in the 33, 41, 61, and 94 GHz 
maps, respectively. The ERCSC catalogue includes 88%, 
63%, 81%, and 95% of the 5a WMAP sources at 30, 44, 70, 
and 100 GHz, respectively. The median of the distribution 
of offsets between WMAP and Planck positions at each fre- 
que ncy are 2.5', 2.1', 1.7', and 1.0' at each frequency (see 
also IPlanck Collaboratior3l2011il for a more detailed discus- 
sion on this subject). Except for the 44 GHz channel, where 
Planck is known to be less sensitive, most WMAP sources 
that failed to be included in the Planck ERCSC (31 sources 
at 30 GHz) are generally at the faint end of the flux density 
distribution (i.e., near the detection threshold) and may 
have flux densities boosted by the Eddington bias or the 
effects of confusion, or may be spurious. The absence from 
the ERCSC of a few brighter WMAP sources (5 sources at 
30 GHz) is probably caused by their variability. 

The Planck-ATCA Co- eval Observations (PACO) 
project fMassardi et al.l [20111 ) has provided measurements 
with the Australia Telescope Compact Array (ATCA) of 
sources potentially detectable by Planck almost simultane- 
ously with Planck observations. 147 ERCSC sources have 
PACO observations at 32.2 and 39.7 GHz within 10 days of 
Planck observations. All these sources are unresolved also 
by ATCA. As illustrated in Fig. [TJ the comparison between 
ATCA and ERCSC flux densities at the nearest frequen- 
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Fig. 2. Euclidean normalized differential number counts at the LFI frequencies. The red circles with Poisson error bars 
show the counts of sources with counterparts in our reference 30 GHz sample. In each panel, the solid curves show the total 
number counts of extragalactic radio sources predic ted by the de Zotti et all ((2005) model. Also shown are: the counts 
estimated at 31 GH z from DASI (grey dashed box; iKovac et al.l ( 2002(1 ) and at 33 GHz from the VS A data (grey box; 
IClearv et al.l (|2005h h t he counts from t he PA CO (grey diamonds; Bonavera et al., in preparation) and t he WMAP 5-yr 
survey s (grey squares; iMassardi et al.l (|2009t )). at the closest frequencies, and the counts estimated by IWaldram et all 
(2007) (grey dashed line), exploiting multi-frequency follow-up observations of the 15 GHz 9C sources. The vertical 
dashed magenta line in the upper panel indicates the flux density completeness limit, 1.0 Jy, estimated for our primary 
sample 



cies (30 and 44 GHz, respectively), shows a reassuringly 
close agreement. The faintest ERCSC flux densities are ob- 
viously enhanced by the effect of the Eddington bias, the 
noise-increased number count of point sources to a given 
detection threshold, which is besides enhanced in the 44 
GHz Planck channel, where the noise level is higher. 

2.3. The 30 GHz extragalactic radio source sample 

2.3.1. Identification of compact Galactic sources 

To minimize the contamination of the sample by Galactic 
sources we have restricted ourselves to \b\ > 5° and we have 
also excluded sources within 5° and 2.5°, respectively, of the 
nominal centres of the Large and Small Magellanic Clouds. 
Outside these regions, a search in the SIMBAD database, 



with a search radius of 16', corresponding to about half the 
FWHM at 30 GHz, has yielded 18 associations of ERCSC 
sources with known Galactic objects (5 PNe, 10 HII regions, 
and 3 SNRs), all within 5' of the ERCSC position. After 
having removed these sources we are left with 533 compact 
extragalactic radio sources detected at 30 GHz, with 97% 
or more of them identified in external catalogues at GHz 
frequencies. This constitutes our primary sample. 



2.3.2. Completeness and uniformity tests 

An indication of the completeness limit of our sample is ob- 
tained by looking at the differential counts (see top panel 
in Fig. [5]) : a sharp decrease of the slope at faint flux densi- 
ties (S v < 1 Jy at 30 GHz) signals the onset of incomplete- 
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ness. Based on the ATCA 20 GHz counts (jMassardi et al.l 
120081: iMurphv et alll2010h . we expect that the slope of the 
counts remains approximately constant over the limited 
flux density range covered by the ERCSC. Therefore, the 
flux density interval (or bin size) containing a fixed num- 
ber of sources must decrease as a power of the central flux 
density. As shown in the top panel of Fig. [31 at 30 GHz this 
happens down to a flux density of about 0.9 Jy, where the 
curve abruptly flattens. 

We need also to test whether the spatial distribution 
of sources in our sample is consistent with being statis- 
tically uniform, as it must be in the case of extragalac- 
tic sources. Deviations from uniformity may be expected 
at lower Galactic latitudes, both because of residual con- 
tamination by unrecognized Galactic sources and through 
the effect of a stronger Eddington bias due to fluctuations 
of diffuse Galactic emission. The bottom panel of Fig. [3] 
shows no significant deviations from a uniform distribution 
on the sky for \b\ > 5° if we adopt a completeness limit of 
~ 1.0 Jy at 30 GHz. Remarkably, the average source density 
at |6| > 5°, D — 24.23 (in sources per sr), is very similar 
to the value found at |6| > 30°, D = 23.71, which guar- 
antees that we are not losing extragalactic sources, in this 
frequency range, when going down to lower Galactic lati- 
tudes, and that the residual contamination due to unrecog- 
nized Galactic sources is negligible. Larger deviations from 
uniformity are found if we adopt fainter completeness lim- 
its. Taking into account both results, we therefore adopt 
S V =1S) Jy as an estimate for the completeness limit at 
30 GHz. Our primary sample is made of 290 sources above 
the adopted S u = 1 . Jy fl ux density limit. As a compari- 
son, iMassardTeTld. (2009) detected 281 sources at \b\ > 5° 
and S v > 1 Jy in their blind survey performed on the 5-year 
WMAP 33 GHz map. 

In addition, we used the Planck HFI frequency channels 
at 143 and 217 GHz to select all the extragalactic sources 
in the ERCSC whose spectra are still dominated by non- 
thermal synchrotron emission at 217 GHz (409 sources with 
a i43 < 0.5; see Fig. 0}. We limit the selection of this sec- 
ondary sample to the above a\ 43 value with the purpose of 
excluding all possible sample contamination coming from a 
second population of sources dominated by thermal dust 
emission (afH > 1.0). The choice of of 43 < 1.0, cor- 
responding to a minimum in the distribution, would not 
change our results since very few sources (< 10) have spec- 
tral indices in the interval 0.5 < 0^43 < 1.0. This selec- 
tion (at above 100 GHz, where there is not yet a complete 
match of ERCSC sources with external catalogs) is useful 
for comparing the outcomes from this (secondary) sample 
with our predictions on the statistical properties of extra- 
galactic sources in our primary sample, selected at 30 GHz. 

3. Number counts 

Figures [2] and [5] show the number counts of extragalactic ra- 
dio sources at the six Planck frequencies from 30 to 217 GHz 
(see also Table [T]). The sharp breaks in the number counts 
at approximately 1.0 Jy (30 GHz), 1.5 Jy (44 GHz), 1.1 Jy 
(70 GHz), 0.9 Jy (100 GHz), 0.5 Jy (143 GHz) and 0.4 Jy 
(217 GHz) signal the onset of incompleteness. 

The results of deeper surveys on small fractions of the 
sky and the WMAP differential number counts are also 
shown, for comparison. The agreement with WMAP data 
is very good. Our differential counts at 30 and 70 GHz are 
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Fig. 3. Completeness and uniformity tests at 30 GHz. Top 
panel: flux density interval, or bin size, including a fixed 
number of sources as a function of flux density. The change 
of slope at about 0.9 Jy signals the onset of incompleteness. 
Bottom panel: uniformity test for sources with flux densi- 
ties S v > 1.0 Jy. The source density, D, in sources per sr, 
within regions at different Galactic latitudes, shows an ac- 
ceptable uniformity. The horizontal grey (solid) and blue 
(dashed) lines show the average source density at Galactic 
latitudes above |6| = 5 and |6| = 30 deg, respectively. As 
it is apparent, the two values are very close to each other, 
within the la level calculated for the overall population in 
the two sky areas considered here, and well inside the la 
normalized Poisson error bars. We also checked that at flux 
densities below about 0.9-1.0 Jy statistically more relevant 
deviations from a uniform distribution begin to appear. 



















IT 







Fig. 4. Spectral index distribution of ERCSC sources be- 
tween 143 and 217GHz. Only sources detected at 143 GHz 
and at 217 GHz have been considered: no upper limits on 
flux densities have been used in this calculation. 
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Table 1. Euclidean- normalized differential number counts per steradian estimated from the Planck ERCSC at 30- 
217 GHz. The 143 and 217 GHz number counts are those represented by blue diamonds in Fig. 5. The bins are centred 
in the log 10 (S„) values and are symmetric in logarithmic scale. 
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somewhat deeper than the WMAP ones at 33 and 61 GHz, 
while at 44 GHz they are somewhat shallower than the ones 
calculated for the 41 GHz WMAP channel. Also our counts 
above the completeness limits appear to join smoothly with 
those from deeper surveys. 

At frequen c ies of up to 100 GHz, the predictions of the 
Ide Zotti et al.l (|2005l ) cosmological evolution model - rely- 
ing on extrapolations from lower frequency data and capa- 
ble of providing a very good fit to almost all data on number 
counts as well as on other statistics of radio sources at fre- 
quencies above 5 GHz - are in generally good agreement 
with our current findings. This result implies that no new 
radio source population shows up at bright flux densities. 
Very few "extreme" or "inverted-spectrum" compact ra- 
dio sources are found in the Planck ERCSC. The emission 
and spectral properties of these sources, which are interest- 
ing in their own right, are di scussed in a companion paper 
(jPlanck CollaborationlfcOllih . 

At higher frequencies (i.e., at 143 and 217 GHz) we also 
plotted the number counts obtained by using the sample 
of radio sources selected at 143 and 217 GHz (blue dia- 
monds; see S I2.3l and Fig.[4]for more details). These number 
counts (calculated from our secondary sample, 284 sources 
at above 0.5 Jy) turn out to be in almost perfect agreement 
with the ones obtained in the previous Section from our pri- 
mary sample selected at 30 GHz (290 sources at S u > 1.0), 
thus confirming that no bright extragalactic radio sources 
are missed by our selection criteria and that the underlying 
parent population turns out t o be (statistica l ly) th e same. 

Figure shows that the Ide Zotti etldl (|2005l) model 
over-predicts the bright counts by a factor of about 2 at 
143 GHz and about 2.6 at 217 GHz, while i t is consis- 
tent with the fainter SPT (jVieira et all I2010H and ACT 
(jMarriaee et al.1 l2010h counts. As discussed in the next 
Section, the discrepancy between the model and our cur- 
rent data is due to a steepening of the spectra of ERCSC 
sources above about 70 GHz, not predicted by the model 
but, at least partially, a lread y suggested by othe r data sets 
(jGonzalez-Nuevo et al.l[200l ISadler et alll2008h . 



An implication of this result is that the contamination 
of the CMB angular power spectrum by extragalactic ra- 
dio sources below the detecti on limit at 143 and 2 17 GHz 
is lower than predicted by the Ide Zotti et all (|2005l ) model. 
Assuming a Poisson distribution (clustering effects are re- 
duced to negligible values by th e very broad luminos - 
ity function o f radio sources, e.g., iToffolatti et~aH ([2005); 
Ide Zotti et al.l (|2010t )) and simply scaling down the model 
counts by the factors mentioned above, the amplitude of the 
angular power spectrum of unresolved sources goes down by 
roughly the same factor 0. This is, however, an upper limit 
to the correction factor, especially at 217 GHz, because if 
we apply the factor calculated above to all flux densities we 
would end u p with a clear under estimate of the faint counts 
measured bv IVieira et al.l (|2010l ). Therefore, if we consider 
much fainter source detection limits, as foreseen for future 
experiments in the sub-mm, the amplitude of the angular 
power spectrum due to unresolved sources stays essentially 
at the same level as predicted by the De Zotti et al. model. 

4. Spectral index distributions 

To study the spectral properties of the extragalactic radio 
sources in the Planck ERCSC we used our reference 30 GHz 
sample above the estimated completeness limit (1.0 Jy; 
§ 12. 3D . Not all of these sources were detected at the > 5a 
level in each of the Planck frequency channels considered. 
Whenever a source was not detected in a given channel we 
replaced its (unknown) flux density by a 5<r upper limit, 
where for a we used the average r.m.s. error estimated at 
each Planck frequency. The upper limits have been redis- 
tributed among the flux density bins by using a Survival 

4 As an example, the amplitude of the power spectrum of radio 
sources, in terms of Ct values, below, e.g., S v < Uy at 217 GHz, 
corresponds to ~ 4.2 x 10 _5 /xK 2 , if calculated from the Planck 
ERCSC differential number counts of Table [T] For comparison, 
undetected radio sources sum up ~ 1.1 x 10~ 4 /xK 2 , if we inte - 
grate the differential number counts of the Ide Zotti et all (2005) 
cosmological evolution model up to S v = 1 Jy. 
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Fig. 5. Euclidean normalized differential number counts at the HFI frequencies (100, 143, and 217 GHz). The red circles 
with Poisson error bars show the counts of sources with counterparts in our reference 30 GHz sample. At 143 and 217 GHz 
the blue diamonds (shifted to the left by half of the bin size, for clarity) show the counts obtained after removing sources 
with 143-217 GHz spectral index indicative of dust emission (see Sect. 2. 3.2). Again, in each p anel, the solid curves 
show the total number counts of e xtragalactic radio sources predicted by th elde Zotti et all (120 05) evolution model. Also 
shown are the SPT (grey squares; IVieira et al.ll2010|) and ACT (grey stars; iMarriage et al.ll2010l counts of radio sources. 
At 100 GHz we also show the estimated counts bv lSadler et al.l ()2008l ) from follow-up observations of a sample of sources 
selected from the 20 GHz ATCA survey (grey dashed line). 



Table 2. Median and standard deviations of the spectral 
index distributions between 30 GHz and the selected fre- 
quency. We adopt the convention S„ oc v a . 



//[GHz] 


44 


70 


100 


143 


217 


median 
error 

a 


-0.06 
0.01 
0.30 


-0.18 
0.01 
0.18 


-0.28 
0.01 
0.17 


-0.39 
0.01 
0.16 


-0.37 
0.01 
0.15 



Analysis technique and, more specifically, by adopting the 
Kapla n-Meyer estimator (a s implemented in the ASURV 
code, ()Lavallev et al.l [1992] )). Since the fraction of upper 
limits is always small (it reaches approximately 30% only 



in our less sensitive channel at 44GHz), the spectral index 
distributions are reliably reconstructed at each frequency. 

Table [5] gives the median spectral indices between 
30 GHz and the other frequencies considered here and the 
dispersions of the spectral index distributions. A moder- 
ate steepening of spectral indices at higher frequencies is 
clearly a pparent. Hints in this directi on were previously 
found by Gonzal ez- Nuevo et al. (l2008|) from their a nalv- 
sis of t he NEWPS sample (jLopez-Caniego et al"1l2007ft and 
also by [Sadler et a l. ( 2008) . Additional evidence of spect ral 
steepening is presented in iPlanck Collaboration! ()2011kf ). 

Our 3 0-100 GHz sp e ctral index is close to the a ~ —0.39 
found bv ISadler et"afl (|2008l ) between 20 and 95 GHz, for 
a sample with 20 GHz flux density S > 150 mJy Moreover, 
our 30-143 GHz median spectr al index is in very good 
agreement with the one found bv IMarriage et ail (|2010t ) for 
their bright (5„ > 50 mJy) 148 GHz-selected sample with 
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Fig. 6. Spectral index distributions for different frequency 
intervals calculated by taking into account all sources se- 
lected at 30 GHz with S v > 1 Jy. There is clear evidence for 
a steepening above 70 GHz. 




Fig. 7. Contour levels of the distribution of a^g 3 vs. o^q 
obtained by Survival Analysis, i.e., taking into account the 
upper limits to flux densities at each frequency. The colour 
scale can be interpreted as the probability of having any 
particular pair of values of the two spectral indices. The 
maximum probability corresponds to air, ~ —0.18 and 



143 
70 



-0.50. 



complete cross-identifications from the Australia Telescope 
20 GHz survey, i.e alf = -0.39 ±0.04. On the other hand, 
(jVieira et all 2010') find that their much fainter synchrotron 
emitting radio sources selected at 150 GHz are consistent 
with a flatter spectral be haviour (mean al 50 ^ —0.1) be- 
tween 5 GHz and 150 GHz. lMassardi et all ((2010ft find mean 
spectral indices af° ~ -0.17 and a\f ~ -0.30 for AT20G 
sources with 150 GHz flux density S > 50 mJy. A flatten- 
ing of the mean/median high-frequency spectral indices at 
flux densities fainter than the ones probed by the Planck 
ERCSC may help to account for the unusually "flat" nor- 
malized counts at 143 and 217 GHz. 

In Fig.[6]we compare the distributions of spectral indices 
over different frequency intervals. There is a clear shift to- 
ward steeper values above 70 GHz: the median values vary 
from agg = -0.18 ±0.01 (a = 0.18) to a^ 3 = -0.52 ±0.01 
(er = 0.22). On the other hand, the distribution of spec- 
tral indices between 143 and 217 GHz is close to the one 
found for a\f (ajU = -0.46 ± 0.01; a= 0.23). This lat- 
ter result is again very similar to the corresponding value 
calculated for all the sources detected at 143 and 217 GHz 
with a-fH < 0.5, i.e., ajYk = —0.51 ± 0.01, as sh o wn in F ig. 
01 Moreover, in the paper IPlanck Collaboration! (|2011kft an 
average value of a = —0.56 (a — 0.29) at HFI frequen- 
cies is found for their sample of 84 bright blazars selected 
at 37 GHz, with flux densities measured in at least 3 HFI 
channels, in full agreement with our present findings. 

Figure[7]presents the contour levels of the distribution of 
ayo 3 vs - a 30 (obtained using Survival Analysis) in the form 
of a 2D probability field: the colour scale can be interpreted 
as the probability of a given pair of spectral indices. 

As already noted, at high Planck frequencies most of 
the extragalactic radiosources are blazars. From the con- 
tour plot of Fig. [7] it is possible to see that the maximum 
probability corresponds to a|g — —0.18 and djQ 3 ~ —0.5. 
A secondary maximum can also be seen at ctjQ 3 ~ —1.2. 
However, a physical interpretation of these features goes 
beyond the purposes of this work and, moreover, more data 



at higher frequencies are needed. A detailed discussion on 
the modelling of the spectra o f this source population is pre - 
sented in a companion paper (jPlanck Collaborationll2011kD . 

5. Conclusions 

We have exploited the ERCSC to estimate the bright counts 
of extragalactic radio sources at 6 frequencies (30, 44, 70, 
100, 143, and 217 GHz) and to investigate the spectral prop- 
erties of sources in a complete sample selected at 30 GHz. 
The counts at 30, 44, and 70 GHz are in good agreement 
with those derived from WMAP data at nearby frequencies. 
The completeness limit of the ERCSC is somewhat deeper 
than that of WMAP at 30 and 70 GHz and somewhat shal- 
lower at 44 GHz. At higher frequencies the ERCSC has 
allowed us to obtain the first estimate of the differential 
number counts at bright flux density levels. At 30, 143 and 
217 GHz, the present counts join smoothly to those from 
deeper surveys over small fra ctions of the sky. 

The Ide Zotti et al.l (|2005ft model is consistent with the 
present counts at frequencies up to 100 GHz, but over- 
predicts the counts at higher frequencies by a factor of 
about 2.0 at 143 GHz and about 2.6 at 217 GHz. This 
implies that the contamination of the CMB power spec- 
trum by radio sources below the 1 Jy detection limit is 
lower than previously estimated. No significant changes are 
found, however, if we consider fainter source detection lim- 
its, i.e., 100 mJy, given the convergence between predicted 
and observed number counts. 

The analysis of the spectral index distribution over dif- 
ferent frequency intervals, within the uniquely broad range 
covered by Planck in the mm and sub-mm domain, has 
highlighted an average steepening of source spectra above 
about 70 GHz. The median values of spectral indices vary 

from a 7 ^ = -°- 18 ± °- 01 ( CT = °- 18 ) to a lf = -°- 52 ± °- 01 
(a = 0.22). This steepenin g accou nts for the discrepancy 
between the Ide Zotti et al.l (|2005ft model predictions and 
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the observed differential number counts at HFI frequencies. 
T he current outcome is a l so in a greement with the findings 
of iPlanck Collaboration! (|201 lkT) on a complete sample of 
blazars selected at 37 GHz. The change detected in the spec- 
tral behaviour of extragalactic radio sources in the Planck 
ERCSC at frequencies above 70-100 GHz can be tentatively 
explained by electron ageing or by the transition to the op- 
tically thin regime, predicted in current models for radio 
emission in blazar sources. However, with present data it is 
not yet possible to clarify the situation. In the near future, 
the data of the Planck Legacy Survey will surely prove very 
useful in settling this open issue. 
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